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Abstract 

We show that the transfer matrix of the A$_ 1 open spin chain with diagonal 
boundary fields has the symmetry U q (SU(l)) x U q (SU(M — I)) x U(l), as well as a 
"duality" symmetry which maps I <-> J\f — I. We exploit these symmetries to compute 
exact boundary S matrices in the regime with q real. 



1 Introduction 



The richness of boundary phenomena is a dominant theme of contemporary theoretical 
physics. Well-known examples include catalysis of baryon decay, D-branes, particle pro- 
duction near event horizons, the Kondo problem, and edge excitations in the fractional 
quantum Hall effect. As much of the boundary phenomena of interest is nonperturbative, 
integrable models with boundaries provide a particularly valuable laboratory for its study 
U, ][J. Magnetic chains associated with affine Lie algebras ||, ||] constitute large classes 
of such integrable models. We focus here on the simplest such class, namely, the Atf_ x spin 
chain. In addition to its value as a simple toy model, this model is of interest in its own 
right: it is related to loop models describing certain self-avoiding walks (see, e.g., and 
references therein), and possibly also to the A$_ l Toda field theory with imaginary coupling 

i, 0. 

Pasquier and Saleur observed || that a certain integrable open XXZ spin chain Hamilto- 
nian has the quantum-algebra (or "quantum group") |J symmetry U q (SU(2)). Kulish and 



Sklyanin later showed [|l(| that this symmetry extends to the full transfer matrix |[L1| . This 



result was then generalized [12] - [14] to higher rank: namely, the transfer matrix of the 
A$_i open spin chain without boundary fields has the symmetry U q (SU(Af)). This model 



was further investigated in Refs. fl5fl, ||16| . Moreover, the most general diagonal boundary 
interactions which preserve integrability were found ||17||,|[L8||. However, the question of what 
symmetry - if any - remains in the presence of such boundary interactions was not explored. 

We show here that by turning on diagonal boundary fields, the U q (SU(Af)) symmetry is 
broken to U q (SU(l)) x U q (SU(Af - I)) x 17(1). Moreover, we find a "duality" symmetry 
which relates / <->• M —I. We exploit these symmetries to compute exact boundary S matrices 
||, which describe the scattering of the model's excitations from the ends of the chain, in 
the regime with q real. For the case M = 2, we recover the results of Refs. [O] and [20 



and for q = 1 we recover the recent results [21 



In Section 2 we review the construction of the open chain transfer matrix, and we exhibit 
its symmetries. In Section 3 we summarize the computation of boundary S matrices. We 
conclude in Section 4 with a brief discussion of some possible generalizations of this work. 
Some technical details have been relegated to the Appendices. 



2 The transfer matrix and its symmetries 

There are two basic building blocks for constructing integrable open spin chains: 
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1. The R matrix, which is a solution of the Yang-Baxter equation 



R12W Riz(\ + A') R 23 (X') = R 23 (X') Ri 3 (X + A') R 12 (X) , (2.1) 

(see, e.g., |2"2"||). We assume that the R matrix has the unitarity property 

R 12 (X)R 21 (-X) = 1, (2.2) 

where i?2i(A) = Vi 2 Ri 2 (X)V\ 2 = Ri 2 (X) tlt2 , t denotes transpose, and Vu is the permu- 
tation matrixP] , and also the property (see also [2B| ) 



with M t = M and 



R 12 (\) tl M 1 R 12 (-\ - 2p) t2 Mf 1 oc 1 , 



[M 1 M 2 , J R 12 (A)] = 0. 



(2.3) 



(2.4) 



2. The matrices K^, which are solutions of the boundary Yang-Baxter equation [23 Q 

R 12 (X 1 - A 2 ) Kt(Xi) #2i(Ai + A 2 ) Kf{\ 2 ) 
= Kf{\ 2 ) R 12 (X ± + A 2 ) Kf(X 1 ) R 21 (X ± - A 2 ) . (2.6) 



The corresponding transfer matrix £(A) for an open chain of N spins is given by fll |, [12 
t(X) =tr M K+(-X-pY To(A) K (X) f (A) , (2.7) 
where tro denotes trace over the "auxiliary space" 0, To (A) is the monodromy matrix 

T (X) = R 0N (X)---R 01 (X), (2.8) 

and T (A) is given by 



T (A) = i?io(A) • • • -Rjvo(A) . 



(2.9) 



■"■The permutation matrix is defined by Vx (8 y = y (g> x for all vectors x and y in an A/"-dimensional 
complex vector space CV- 

2 We use here a K + which is related to the K^ ld used in earlier work Q - Q by 

K+ d (X) = MK+(-X - pY . (2.5) 

The matrix K^ ld satisfies an equation more complicated than ( |2.6| ) , but the expression for the transfer matrix 
in terms of K^ [d is simple. The matrix K + which we use here obeys the same boundary Yang-Baxter equation 



as K , but the expression for the transfer matrix in terms of this K + (see Eq. (2.7)) is more complicated. 
The two approaches are of course equivalent. 
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(As is customary, we usually suppress the "quantum-space" subscripts 1 , . . . , N.) Indeed, 
it can be shown that this transfer matrix has the commutativity property 



[f(A),t(A0] = O. 



(2-10) 



In this paper, we consider the case of the A$_ x R matrix [25 
Ri'iXjjj.jj = 1, 



R\2\X)jk,jk 
Rl2\X)jk,kj 



sinh(— irjX) 
sinh (r}(— iX + 1)) 

sinh?7 
sinh (r](—i\ + 1)) 

l<J,k<M, 



JriXsign(J-k) 



(2.11) 



which depends on the anisotropy parameter rj > 0, and which becomes SU(J\f) invariant for 
77 — > 0. This R matrix has the properties ( |2.2j ) and 



, with |L7 

M jk = 5 jk e^- 2l+l \ p = iM/2. 
Moreover, we consider the N x M diagonal K matrices given by []TB] 

K-(X) = tf (0 (A,£_), 
K + (X) = K {l) (X,£ + -^), 



(2.12) 



(2.13) 



where 



#(q(A,0 = diagi^ 




a = sinh (77 (£ + z'A)) e 



b = sinh (r/(£ — z'A)) e 



(2.14) 



for arbitrary £, and any Z 6 {1 , . . . , JV— 1}. Eq. (|2.14 ) is the most general diagonal solution 
of the boundary Yang-Baxter equation with the A$_ 1 R matrix, up to an irrelevant overall 
factor. 

We shall usually denote by £(z)(A , £_ , £+) the corresponding open spin chain transfer 
matrix Q 

t (0 (A , e- , £+) = tr M K (0 (_A - p , e+ - y) ^(A) tf (0 (A , £-) To (A) . (2.15) 



3 The more general transfer matrix i(; + ,;_)(A,£_ ,£+) constructed with -f^(; T )(A,£) also forms a one- 
parameter commutative family. For simplicity, we consider here (as in [pi]) the special case l + = I- = I. 
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The corresponding open spin chain Hamiltonian Ti. open is related to the derivative of the 
transfer matrix at A = 0: 



open 



4 sinh(^-) tr MK {l) (-p , £+ - f ) d\ 



A=0 



i d 

AT-l 



logtr MK 



(IV 



d 



2 ' 



A=0 



5 Hnn+1 + 4sinh(r / £_)Z\^ ) l(A ' ^ 

tr M K (l) (-p^ + -f)n m 
trMif (/) (-p,£ + -f) 



where the two-site Hamiltonian 7^ is given by 



A=0 



A=0 



One can verify that the Hamiltonian is Hermitian. 

The parameters £~ x may be regarded as certain boundary fields. For £ 
matrix Q2.14j ) evidently becomes proportional to the identity matrix 

WA,£)~ \e* 1, 



(2.16) 



(2.17) 



oo, the K 



(2.18) 



and so the transfer matrix is U q (SU(Af)) invariant [13], |14]]. We now exhibit exact symme- 
tries of this transfer matrix for finite values of £ T . 



2.1 Duality 

The duality symmetry of the transfer matrix is a remnant of the cyclic (Zj^) symmetry [^5 
of the A$_i R matrix. Indeed, the R matrix ( |2.11|) satisfies 



Ui R 12 (X) ur 1 = V 2 (-\)- 1 R 12 {\) V 2 {-\) 
U 2 R 12 (X) U 2 l = V^X)- 1 R 12 {\) V x {\) , 

where U is the J\f X J\f matrix 

Ujk = <5j,ai-i + <Wj'^i,k 
which has the property = 1, and V(X) is the matrix 

V(\) jk = e- ir > x 6 jtk - 1 + e ir > x 6jtfj6 1>k . 



(2.19) 



(2.20) 



(2.21) 
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The corresponding quantum-space operator U defined by 

U = U l U 2 ---U N (2.22) 
therefore has the following action on the monodromy matrices 

UT {\)W l = Vo(X)- 1 T (X) V (X) , 

U f (A) U- 1 = Vo(-A)- 1 fb(A) V (-X) , (2.23) 



and the transfer matrix ( 2.15 ) transforms as follows: 

W l f (0 (A,£_,£ + )Zr' = tr {^(-A)'Mir (0 (-A-p,e + -y) V(A)- ; ) 

x r (A)(y(A) , if(o(A,e-) V(-A)- ; )oT (A)}. (2.24) 



We now observe that the K matrix ( p. 14 ) satisfies 



V(\) 1 K (l) (\,t_)V(-\)- 1 = -e-^-^cv-oCA.-e-), 
V(-A)' MK (0 (-A-p,e + -y) V(A)- ; = _ e ^-iM e #-% % _ j)( _ A _ /)i _ (+ + ^ ) 

Z = 1,... ,N~ !■ (2-25) 
We conclude that the transfer matrix has the "duality" transformation property [] 

U l t {l) ( A , £- , e+) = e"^-^ t (10 (A , ^ , , (2.26) 

where 

= -<£+ + a/" , 

Z' = M -l . (2.27) 

Notice that U l 'lA l = 1. The transfer matrix is "self-dual" for £_ = = 0, £+ = C+ = A/"/2, 
and 1 = 1' = M/2. 

2.2 Quantum algebra symmetry 

In the defining representation of SU(Af), we identify (following the notations of our previous 



paper pl| ) the raising and lowering operators 

3 +(k) = e k , k+ i , r {k) = e k+1 , k , k = \,...,N-\, (2.28) 

4 The prefactor e''^ -2 ^ can be absorbed in the definition of K + , and is not significant. 
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which correspond to the simple roots, and the Cartan generators 



k = 1 



,7V- 1 



(2.29) 



where are elementary TV x TV matrices with matrix elements (ek,i) ab = Sk, a 8i,b- We denote 
by , s^> the generators at site n, e.g., 



e ] = i 



n. th 
I 



1® s (k) ®1 ® 



n = 1,... ,iV. 



(2.30) 



The corresponding generators of the quantum algebra U q (SU(Af)) which act on the full 
space of states are given by (see, e.g., |§,[[jj| ) 



N 



j±(k) = £ 



(*0 



s w 



n=l 
N 

£ 

n=l 



,(*) 



fc = 1,... ,7V- 1. 



7 4 fc) /2 



(2.31) 



These generators obey the commutation relations 



J+{k) J-(j) 



5 k 



k,j 



s (k) 



S (k) ) J+ 0-)l = (24j _ ^ _ 5k+i j) j+U) 5 (2.32) 



where [x] = (q x — q x )/(q — q 1 ). 

We claim that the transfer matrix t(j) (A , £_ , £+) has the invariance £/ g (.S'Z/ (/)) x U q (SU (TV 
£7(1). In particular, we shall now show that 



t (0 (A, £_,£+),£« 
t (0 ( A, £_,£+), 



0, k = 1,... ,7V- 1 
0, fc^I, 



(2.33) 
(2.34) 



where q = e~ v . Readers who are interested primarily in the calculation of boundary S 
matrices may now wish to skip directly to Section 3, which may be read independently of 
the proof of Eqs. ( |2.33| ) and ( |2.34| ) that is presented below. 

It is convenient to introduce, following Sklyanin pT| , the quantity Tm (A , £) defined by 

T {1) (A , = T (A) K {1) (A , T (A) , (2.35) 

in terms of which the expression ( |2.15| ) for the transfer matrix acquires the more compact 
form 



TV 

t (0 (A , £- , C+) = tr M K (l) (-A - p , £+ - y) 7^ (A , £_) . 



(2.36) 
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The proof of the first relation ( |2.33|) is straightforward. From the identity 



,(*) 

'2 j 



it follows that 

y k \T (x) 

and similarly for To (A). Hence, 
where we have also used the fact that 



4 fe) ,%o(A,o 



(2.37) 



(2.38) 



(2.39) 



tr n 



0. We conclude that 



S« , M Jfy (-A - p , C+ - — ) T (l) (A , £_) 



tr 



4 fc) , M # (0 (-a - p , e+ - y) ^ } (a , e_) 



AT, 



0. 



(2.40) 



where the final equality follows from the cyclic property of the trace. 

Our proof of the second relation ( |2.34| ) relies on the well-known fact that in the limits 
A — > ±?oo, the monodromy matrix becomes an upper/lower triangular matrix whose elements 
can be expressed in terms of the U q (SU(J\f)) generators. (See, e.g., 0,|13[, Hi~6| .) It is 
therefore convenient to introduce notations for the corresponding limits of the R and T 
matrices: 



for A — > ±zoo, with 



Schematically, 



J- n 



#i 2 (A) 
To(X) 

r (A) 

rp± _ 

1 o — 
f± - 

( * J- (1) 
* 





p=RD± 



± 

■0 J 



Ron ' ' ' R01 > 

-^10 ■ ■ ■ R%o ■ 



(2.41) 



(2.42) 



* 



T, 



> 



7 



To ~ 







i 



(2.43) 



where the diagonal matrix elements involve only the Cartan generators. We shall not need 
more explicit expressions for Tq and T ± , which can be found e.g. in ||16| . 

It is also important to realize that in these limits the K matrices ( [2.14| ) become projectors: 



for A — > ±zoo, where ILa are the orthogonal projection operators 



(2.44) 



n 



(0 



diagil 



,1,0, 



M- 



n 



(i) 



diag[0 ,0,1,... ,1 



(2.45) 



Af-l 



and therefore Tit*. + IL* = 1. 



It follows that for A — > ±zoo, the quantity 77n o(^>0 defined in Eq. ( |2.35| ) tends (up to 



an irrelevant factor) to T^; , where 



2 (i) o — J o ii (o o J o • 



It is easy to show (see Eq. ( |A.4| ) in Appendix A) that 7^t and 7^ obey 



-(/> 



n 1 * 1 T ± n 1 * 1 - T 1 * 1 
ii (Z) 2 (l) o o — 



(2.46) 



(2.47) 



With the help of Eq. (|2.43| ), one can also see that T^X depends on the generators J ±( - k ^ 
with k < I, while depends on the generators J ±<yk ^ with k > I. 

To prove the relation ( [2.34 ), it suffices to show that the transfer matrix commutes with 



'(0 



where are defined by 



*(i)(A,f-,£f) ,r, 



(2.48) 



(2.49) 



and P 1 * 1 are M x jV matrices which obey 



n « p± nfo 



(2.50) 



but which are otherwise arbitrary. Indeed, with suitable choices of P + and P~, one can 
project from and 7^ the generators J ±f - k ' with k I. For example, with P + 
one projects out J ±( - k ^ with k < I. 



;±(&) 



The quantities TVn evidently have a structure similar to that of the transfer matrix. The 



proof of Eq. (|2.48|) is similar to the proof of the fundamental commutativity property fl2.10|) 
of the transfer matrix. Indeed, following Sklyanin, one can show using the properties 



( pT4l) that 



r±t (0 (A,e-,e + ) = t ri (p^j j* 1 tr 2 (M 2 /f (02 (-A-p,e + -y)% 2 (A,e-)) 



Similarly, 



t (0 (A, £_,£+) 7"(f) 



tr 12 {i^ P ± ^ M ^R^K {1) 2 (-A - p , e + - 



x 



(2.51) 



tr 2 (m 2 /% 2 (-A -p,t+- y)2( 2 (A , £_)) tn (Pf^ J* 1 



tri 



12 1 ^(/) 2l 



Af 



x 



(2.52) 



It is easy to see that the quantities in square brackets in Eqs. ( [2.51D and (|2.52|) are equal. 
Indeed, since 7(z)(A,£) obeys the boundary Yang-Baxter equation 



Pi2(Ai — A 2 ) i(Ai , £) P 2 i(Ai + A 2 ) T(i) 2 (A 2 , £) 
= 2(A2 , Pi2(Ai + A 2 ) Ty) x(Ai , £) P 2 i(Ai — A 2 ) , 
by taking Ai — » ±zoo we obtain the relations 

R 12 ^Z) i R fi T{i) 2 (A , = 7(J) 2 (A , 0-^12 1 -^21 • 



(2.53) 



(2.54) 



In order to prove the commutativity (|2.48|) , it therefore suffices to show that 

tr (AC) = tr (PC) , 

where 



A = Rf 2 P^M^R^K (l)2 (-X-p^ + - 



(2.55) 



B = K (l)2 (-X-p^ + -^)Rf 2 Pt tl M^ 1 Rl, 



C 



R?2 % i T Q) a(A , £-)MiM 2 . (2.56) 



tl2 



We establish this result by showing that A — B and C are trace orthogonal matrices. That 
is, we show that for some projection operator n, 

ucu = c 

n {A-B) n = 0, (2.57) 
from which it immediately follows that 

tr [(A -B)C}= tr [(A - B) II C II] = tr [II (A - B) II C] = . (2.58) 
Indeed, in Appendix A we show that 

nf l)1 (Rf 2 T± 1 Ri 1 T {l)2 (\,OM 1 M 2 )llf l)1 = Rt 2 T { t )1 Rf i r {l)2 (X,0MiM 2 , (2.59) 

nfo i {R^p^m^r^k® 2 (a , o) nf x = nf 1 {k (1) 2 (a , or^p? 1i M{ 1 r^) nf x , 

(2.60) 

thereby completing the proof. 

We have considered here symmetry transformations which leave the anisotropy parameter 
7] unchanged. In Appendix C we briefly discuss related transformations, which however 
transform 77 — > —77. 



3 Boundary S Matrix 

The eigenstates of the transfer matrix tm{\ , £_ , £+) can be constructed by the so-called 
nested Bethe Ansatz using the pseudovacuum un^ defined by 

= V(fc) ® • • • ® , fc = 1 , . . . ,JV, (3.1) 

V v ' 

A' 

where t> a.) are the standard A/"-dimensional Cartesian basis vectors with components vn^ a = 
5k, a- The corresponding Bethe Ansatz equations (BAE) can be cast in the formf] 

frjg(Ag>) = jy, a = l,-,J^, J = l,-.. ,/V-l, (3.2) 

We restrict to real Aq , which suffices for computing boundary S matrices. 
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where Tim (A) is the so-called counting function (explicit expressions are given below), and 
{^a } are i n t e g ers lying in certain ranges which serve as "quantum numbers" of the Bethe 
Ansatz (BA) states. We label the holes in the ranges by {J^} , a = 1,... ,u^\ The 
corresponding hole rapidities {A^P} are defined by 

= J®, a = l,...,i/W. (3.3) 

The nature of the ground state and excitations is similar to that of the corresponding 
isotropic model |27 |, [p8l,|21|, except that Lie algebra symmetries are now replaced by their 



g-deformations. Indeed, the ground state is the BA state with no holes, i.e., with M — 1 
filled Fermi seas. Moreover, the BA state with one hole in the j th sea (i.e., v^°' = 5k j 
for k — 1 , . . . , J\f — 1) is a particle-like excited state which belongs to the fundamental 
representation [j] of U q (SU (AT)), corresponding to a Young tableau with a single column of 
j boxes, as shown in Figure [1]. This excitation carries energy s^(A^) and "momentum" 



> 3 



Figure 1: Young tableau with a single column of j boxes, corresponding to 
the fundamental representation [j] of U q (SU(J\f)) 



pU'(\V'), where 

d 



— p U) (\) = 2tts (j) (A), (3.4) 
dX 



and the Fourier transform of s^(A) is given byQ 



{j) _ sinh((AT-j)r/|fc|) 
S W ~ sinh(A^|) • [6 - b) 



6 Since the boundary interactions break the bulk U q (SU(J\f)) symmetry, states should - strictly speaking 
- be classified according to the unbroken symmetry U q (SU(l)) x U q (SU(Af — I)) x £7(1). However, one 
expects that the effects of the boundary should be "small" at points of the chain that are far from the 
boundary. In particular, bulk multiparticle states should "approximately" form irreducible representations 
of U q (SU(M)), and therefore we refer to such states by U q (SU(J\f)) quantum numbers. 

7 Our conventions for Fourier transforms are as follows: 

/(A) = !L E e- 2 ^ kx f(k) , f(k) = / dX e 2 ^ kx f{\) . (3.5) 
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Introducing the lattice spacing a, and taking the continuum (a — > 0) and isotropic (77 — > 0) 
limits with -e ^» = // finite, one obtains a relativistic dispersion relation, with mass 
m^) = /i sin (jrj/J\f). 

We define the boundary S matrices r for a particle of type [j] by the quantization 
condition |3(J 

(^(^5+ w ^ w _i)|F)) = 0. (3.7) 

We compute these S matrices from the ■h terms in the distribution of roots of the BAE for 
the corresponding BA states |A^}. Indeed, we rely on the key identity 

;>(A) + ^(A)-2 S « (A) .^» (A) , (3.8) 
where the "density" cr$(A) defined by 

< w = JiTx h ^ w (3J) 

describes the distribution of roots of the BAE. This method |Kj is a generalization of the 
approach pioneered by Korepin |3l] and Andrei-Destri |3l| for calculating bulk two-particle 
S matrices. (See also |33| .) 

In order to compute each of the boundary S matrix elements, we shall need the densities 
o"m (A) corresponding to each of the pseudovacua up) , . . . ,00^. It can be shown that the 
Bethe Ansatz equations (and hence, the counting functions fo[n (A) and densities 0m (A)) 
corresponding to the first / pseudovacua up) , . . . , uq) are all the same; and the BAE corre- 
sponding to the last Af — I pseudovacua uy + i) , • • • , U(j\f) are ah the same. This result is quite 
plausible, in view of the U q (SU(l)) x U q (SU(Af — I)) symmetry of the transfer matrix. For 
completeness, we provide in Appendix B a proof of this result. 

Corresponding to the first I pseudovacua wm , . . . , wm, the BAE for the transfer matrix 
*(J)(A , £- > 6+) are given by |Tgj 



e 2 ^_ +z (A«) e_ (2C+ _ 2 ^ +0 (A«) + (1 - ^ 

x n e -i(^ } - A r i} ) e -i( A ? + x i 1] ) n ^ - a ?) ^ + 

/3=1 /3=1 

(3+ a 

x II MAW - Ag +1) ) + Ag +1) ) 

0=1 

a = l,...,M^, j = l,... ,JV-1, (3.10) 
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where 



e n (A) 



sinr7(A + f) 
sin?7 (A — y) 



(3.11) 



and M (0) = JV , M (A ° = , X { ° ] = A^° = . The requirement that solutions of the BAE 
correspond to independent BA states leads to the restriction < A^ < £. For simplicity, 
we restrict to "weak" boundary fields > \{M — 1), > Af — \. The energy as a function 
of and £ + is discussed for the case M = 2 in 



Defining the counting function so that the Bethe Ansatz equations take the form (J3] 
we find 



271 

Af(j-i) 



+ E [^(A-A^ + g^A + A 



-g2£_+z(A) + g 2 $ + -2^+z(A)J 



p=i 



E [g 2 (A-A^)+g 2 (A + A^) 
0=1 



+ E [?i(A-A? +1) )+?i(A + A« +1) )]} 

/3=1 



(3.12) 



where <7 n (A) and r n (A) are odd monotonic-increasing functions defined (for n > 0) by 
q n (X) = n + i\oge n (X) , 

7T 7T 7T 

-Ti < q n (X) < 7i for - — <A<— , and q n (X + -) = q n (X) + 2n , (3.13) 

IT] IT] T] 



r n (A) = ilog£ n (A), g n (X) 



COST] 



COST] 



(A-f)' 



-TT < r n (A) < 



7T for - — <A<— , and r„(A + -) = r n (A) + 2tt . (3.14) 
2?7 2?7 7] 



We consider now a multi-hole state, with the number of holes in each of the seas given by 
, i/^ -1 ). In the thermodynamic limit, we obtain a system of linear integral equations 



for the densities ( |3.9|) by approximating the sums in fain (A) by integrals using p0 



^ E ^ } ) 

iv a=l 

In this way we obtain 



ir/2ri 



g(X') a$(X') dA'-lE^(Ag))-^ 



1 

2N 
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J(0) +J (^) 



(3.15) 



m=l 



E ((* + * *fo°0 ( A ) = 2ai(A)^! + -{a 2 (A) + 6 2 (A) 

13 



+ (di(A) + 61(A)) (-1 + 5 jt i + + (~ a 2^+l(X) + a 2 ^ + -2Af+lW) 

+ EE(^-tt m + /C(A + A( m) ) Jm )}, j = (3.16) 



m=l a=l 



where cr[n S (A) is the symmetric density defined by 



4lM 



* denotes the convolution 



4? (A) A>0 
($(-A) A<0 



tt/27? 



(3.17) 



and 



rn 

a„(A) 
&n(A) 



(J*9)W= f(X-X')g(X')dX' 

-n/2r] 



o-2{X)S mt j — ai ( A) (<5 TO j_i + 5 mj+ i) , 
1 dq n (X) r] sinh(?7n) 



(3.18) 



2tt dX tt cosh(^n) — cos(2?7A) ' 

1 dr n (X) r] sinh(?7n) 
27r dA 7r cosh(^n) + cos(2r^A) " 2i] 



71 



(3.19) 



Using Fourier transforms, we obtain the solution 



a 



(0 • 



(A) = 2s<fi(\) + 6<T$(\) 



(0 



E (^jm * [02 + 62 + (01 + 61) (-1 + <W + <W_ X )]) (A) 



Af-1 i/( m ) 



+ E E - Ai m) )^> - fc iTO (A - AL m) ) + (A^ m) - -Ai m) )] } , (3.20) 



m=l a=l 



where lZj m (X) has the Fourier transform (see, e.g., pq| ) 



e ??]*| gjn^ (j77. < r]|fc|) sinh ((A/" — m>)?7|fc|) 



sinh (Afi]\k\) sinh (^|/c|) 
with m> = max(m , m') and m < = min(m , m'), and 

5a (i)i X ) = * [~ a 2i-+i + a 2i+ -2M+i)) (A) 

has the dependence on the boundary parameters £ T . We note also 

a n (k) = , b n (k) = (-) k a n (k) , n > 



(3.21) 



(3.22) 



(3.23) 
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The above density a^) S (A) corresponds to the first I pseudovacua cj(i) , . . . , LUyy We also 
need the density cr^ s (A) corresponding to the last J\f — I pseudovacua (jl>(i+i) , . . . ,0J(M), as 
already remarked. We recall that under duality, the transfer matrix transforms as ( p. 26 ) 

U l t (0 ( A , i_ , £+) U~ l oc t (n (A , (3-24) 

while the pseudovacuum transforms as 

U l U(k) = U(k-i) , (3.25) 



where we identify U(k) = u^+M)- It follows JH]] that the BAE for £(j)(A , £_ , £+) with pseu- 
dovacua uj(i + i) , • • • , wr/v) are the same as the BAE for tn')(X , , £+) with wm , . . . , wm; i.e., 
the BAE are given by Eq. ( |3.10|) , except with £ T — > £' and / — > I'. It follows that the 



corresponding densities <7^ S (A) are given by Eq. ( |3.2U| ), except with 



6(T i)\ x ) = { n i,N-i * (a 2 £_-M+i - a2( + -M+i)) (A) • (3.26) 

For simplicity, we compute boundary S matrices only for the cases [1] and [Af — 1] = [I], 
for which the S matrices act in the Af- dimensional complex vector spaces CV and CV, 
respectively. 

We begin with the case [1] , for which the boundary S matrices ^ are diagonal AfxAf 
matrices of the form 

s (i) [i] = diagia^ , . . . , af t) , pfa , . . . , (3^) . (3.27) 

S v ' v ' 

I N-l 

Choosing the state (A^) in Eq. (3/7) to be the BA state constructed with any of the first 
I pseudovacua , . . . ,wm having one hole in sea 1 (i.e., = t^i), it follows from the 
identity ( |3.8| ) that (up to a rapidity-independent phase factor) 



"(0 a (i) ~ 



exp I^ttN £ } (a$(X) - 2s (1) (A)) dx\ , (3.28) 



where crf z V(A) is given by Eq. ( |3.20|) with j = 1. We now observe that 



£ { ) [1I (A - A (1) ) + TZ (A + A (1) )] d\ = jT ' 2K (2A) dX . (3.29) 
Moreover, we note the identity 

^ i-e- 2 ^ A; r (?2 ( / u + ( 5)r g 2(/i + 7) 
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as well as the q- analogue of the duplication formula [35] 



[l + q) 



2x-l 



r„2(x) Y r ,2(x + -) 



(3.31) 



where T q (x) is the g-analogue of the Euler gamma function. 
We find f\ 

q -i2~xma+i)/« So{ xW) 



a 



x 



T q 2Af 



(ir (e- + 1 - 1 + iA«)) r> (e_ + M - \ - iAW)) 



+ z - i - iAd: 



(0 



g -«A«(l-i+2AO/^ o( X(l)) 

r^(^(e + -^+z-| 



r> (£ (e-+^-| + zA(i))) 



aw 



-§ + iA« 



where the prefactor So (A) is given by 

r> - 1) - r> (at + *a 



| - zA« 



(3.32) 



5 (A) 



U (K^ " 1) + ^ I> - zA 



(3.33) 



Furthermore, choosing the state |A^) in Eq. ( |3.7| ) to be the BA state constructed with 
any of the last Af — I pseudovacua , . . . , W(aT) having one hole in sea 1, we obtain the 

relation 



a (i) a (i) 



exp UnN £ W (o§\\) - <$(A)) d\\ . 



(3.34) 



We note that 



<$'(A)-<$(A) = ^(A)-^W(A) 



'(0 



(0 



'(0 



(a 2 g__i(A) - a 2 ^ + -2N +2«-i(A) 



(3.35) 



where and are given by Eqs. ( |3.22| ) and ( |3.26| ), respectively. We conclude 



P(0 



+ 

(i) 



-e2|+-2AT+2Z-l (A (1) 



(3.36) 



8 Evidently, Eq. ( 3.28| ) determines the product OL%.otny Further conditions are necessary to determine 
ajjs and a%, separately. Following |p0| , we assume that the parts of and afn which do not originate 
from 8a^ are equal. 
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where we have resolved the sign ambiguity by demanding that the S matrix be proportional 
to the unit matrix for A'- 1 -' = 0. 

We remark that there are relations among the boundary S matrices. Indeed, the above 
results for S^s m are valid in the regime > |(jV- 1), > M- i, to which we now refer 
as regime I. Corresponding results for regime II defined as < — |(jV— 1), £+ < — Af + \ 
can be obtained by a similar analysis, starting with the BAE ( |3.1U| ) and replacing — > 
— , £ + — > — £ + + A/". We observe that the S matrices in these two regimes are related by 
duality: 



a 



t (J) 
(0 



P{i) 



t (//) 
(0 

T (J/) 

(0 



(3.37) 



where £'± and /' are given by Eq. (|2.27| ). Moreover, there is a relation between Sn\ m and 
^(O [i]' 



D (0 [i] 



q+ {II) 



(3.38) 



Finally, we consider the case [Af — 1]. The boundary 5 matrices S 1 ^ ^ ^ are diagonal 
M x A/" matrices of the form 



" v ' V v ' 

/ N-l 

For this case we must consider one hole in sea M — 1. Noting that 
At) % ~ exp j^TTAT j( «- 1} (A) - 2s^- 1 \X) 



(3.39) 



(3.40) 



we obtain 

a- _ ^-i2A<^- 1 )(l-i+A r 

P(0 - 1 



T„2Af 



x 



+ 



a+ _ „-i2A( AA - 1 )(l+i+AT 
P(0 _ V 



X 



r„2A^ 



+ |(A~ - 1) - g^)) I> (g_ + |(A~ - 1 ) + zA^- 1 )) 



+ i(A" - 1) + *A(^-D)) i> (i (e_ + ±(Af - 1) 
- |(A-+ 1) I> (g + - |(A-+ 1) - 



-iA^- 1 ))) 



iA(^-i) 



- i(A/-+ 1) - tAC^-D)) T q ^ (T ( e+ _ i (A r+ i) +,A(^-i))) 



(3.41) 
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where So (A) is given in Eq. ( j3.33|) . Moreover, 



a 



(0 
+ 
SL 

(0 



-e 2e _ +2 z- A r-i(A (Ar - 1) ) 
-e^-Ar-UA^ 1 )). 



(3.42) 



For the case M = 2, we recover the results of Refs. and f20 |. In the isotropic limit 
i] — > 0, we see that g — > 1 and r 9 (x) — > r(x); hence, we recover the recent results [2I|. 



4 Discussion 

The transfer matrix t (l) (A,£_ ,£+) has the symmetry C/ g (S£/(Z)) x £/, (SU(M-l)) x 17(1) 
as well as a duality symmetry I <-> A/" — /, as we have shown in Eqs. (|2.33|) , fl2.34|) and 
( p.26| ), respectively. It should be possible to show, generalizing |13[, (T^j, [j3"B| , that the 
corresponding Bethe Ansatz states | ) are highest weights of f/ g (SU(l)) x C/ g (SU(J\f — I)). 
That is, 

j+( fc )|) = o, fe^r (4.i) 

This symmetry was used here for computing boundary S matrices. We expect that this 
symmetry will also be valuable for computing other quantities of physical interest (e.g, 
form factors, correlation functions, etc.) and for elucidating the phase structure of the 
corresponding vertex model. Although we have considered here the noncritical regime, we 
expect that this approach can be extended to the critical regime |g| = 1. 

The "mixed" boundary condition case l + ^ /_ merits further investigation. One expects 
that the boundary S matrix for one end of the chain should be independent of the boundary 
conditions at the other end. However, for mixed boundary conditions, the unbroken sym- 
metry group of the transfer matrix is smaller, and therefore, the arguments presented here 
require further refinement. 

We have seen that by turning on diagonal boundary interactions, the symmetry algebra 
of the transfer matrix is broken to a nontrivial subalgebra. We expect that this is a general 
phenomenon. More explicitly, consider the R matrix associated with the affine Lie algebra 
9 [01) 0- It is known [13] that the corresponding open spin chain transfer matrix with 



K — 1 (i.e., without additional boundary interactions) has the symmetry U q (go), where go 
is the maximal finite-dimensional subalgebra of g . We expect that the transfer matrix with 
diagonal K ^ 1 will have as its symmetry a nontrivial subalgebra of U q (go). 



18 



We emphasize that we have considered here only diagonal boundary interactions. It would 
be interesting to look for symmetries of the transfer matrix with nondiagonal K matrices. 
Unfortunately, even for the simplest case the Bethe Ansatz solution is not yet known, 
which precludes computation of the boundary S matrix in the manner described here. 

We also emphasize that here we have restricted our attention to exact symmetries of 
the finite chain. For the (semi) infinite chain, one expects additional affine symmetries 
|19]|,||37||. In field theory, such symmetries are generated by fractional-spin integrals of motion 
33|,[3£],0, which can exist even for nondiagonal boundary interactions [40 . 
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A Proof of Eqs. ( ggg ) and ( gggj ) 



We prove here the relations ( [2.591) , (|2.6Q|) which enter our proof of quantum-algebra symmetry 
of the transfer matrix. 



i Kyi n«) i — R\2 i , (A.i, 



We begin by noting the identities 

'■(O 1 -^12 1 — 1L 12 1 J 

which imply 

np) i Rfi x = l r 21 , (A. 2) 

and therefore 

TT± rp± jt± rp± -pr± 

xi (Z) 1 xi (Z) — J U (Z) ' 

H(Q ^0 = H(t) ^0 • (A-3) 

It is now easy to show that 

L (z) z (z) ~~ • 

Indeed, recalling Eq. ( |2.46|) , 



%V\ — %X\ ' (A- 4 ) 



Hp) ^(o ~~ o Tq t — f — 7^ . (A. 5) 
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In a similar manner one can prove that 

n (0 i (^12 7(0 i R fi 7(0 2(A , 0) x = Rf 2 T^j x R 21 7(0 2(A , 
This relation, together with the fact 

^%MUf l) = Uf l) M = MUf l) , 



imply the first desired result ( |2.59|) . 
We next note the identities 



which imply 



H(o i R12 i — x Rf 2 , 



n^) i R21 ! — R 2 \ 1 . 



(0 



Let us denote by R' T the "projected" R matrix, i.e., 

= ! Rf 2 ! . 

It is easy to see that this projected R matrix commutes with Kq) 2 (A , £), 
Indeed, we see from Eq. (|2.14 ) that 

^ {0 (A,o = anJ ) + 6n w . 



(A.6) 



(A.7) 



(A. 



(A.9) 



(A.10) 



(A.ll) 



(A.12) 



Since 



-^12 5 



± 

(0 i Xi (0 2 



0, then 



^12 > 2 



0, and the result (|A.11| ) for R' 12 immediately 



follows. The proof for R' x \ is obtained in similar fashion. 
All the pieces are now in place for proving the relation 

n± x (rtp^m^rtk^ 2 (a , 0) n± , = (k {1) 2 (a , qb&p? ^m^r^) n± t , 

(A.13) 

which is the second desired result ( |2.60| ). Indeed, 

LHS = R^Pt^M^RlK^ 2 (A , = #(Q 2 (A , ^P^Mf 1 ^ = RHS . (A.14) 
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B Bethe Ansatz Equations for general pseudovacua 



A priori, one expects that the Bethe Ansatz equations (BAE) depend on the choice of pseu- 
dovacuum. We show here that the BAE for the transfer matrix t(/)(A , £_ , £ + ) corresponding 
to the pseudovacua um , . . . , um are all the same; and that the BAE corresponding to the 
pseudovacua , ■ ■ ■ , &(Af) are all the same. 

We first observe that the operators ( J*^) , where N is the number of sites, transform 
one pseudovacuum into another: 

(j +ik) Y w w = ^,4^1. (B.l) 

The U q (SU(l)) x U q (SU (Af — I)) symmetry of the transfer matrix ( |2.34|) implies the (weaker) 
result 



t (0 (A,e-,e + ),(^ ±(fc) ) 



N' 



0, k^l. (B.2) 



Let A/» | (A , £_ , ) denote the eigenvalue of the transfer matrix corresponding to the first 
pseudovacuum 

(a , e_ , e+) w ( i) = (a , e- , e+) ^ ■ (b.3) 

Acting repeatedly on both sides of this equation with the lowering operator (j~( k ^ , we 
conclude that the transfer matrix has the same eigenvalue with each of the first I pseudovacua, 

%) (A , £_ , £+) cu (k) = A$ (A , £_ , £+) w (fc) , fc = 1 , . . . , I. (B.4) 

Similarly, let A'^(A,£_ , £+) be the eigenvalue of the transfer matrix corresponding to the 
last pseudovacuum coy^, 

%)(A , C- , 6+) ^(AO = AjJ^A , £_ , 6+) . (B.5) 

Acting on both sides with the raising operator (j + ( k A , we see that the transfer matrix has 



the same eigenvalue with each of the last M — I pseudovacua, 

t (0 (A , , £+) oj {k) = A' ( ( ° } (A , , e+) oj {k) , fc = / + l,...,A/\ (B.6) 

We have shown so far that the vacuum eigenvalues of £(z)(A , £_ , £+) with u;^) , . . . ,u>^ 
are the same, and those with 0J{i+\) , . . . , are the same. The eigenvalues of £(j)(A , £_ , £+) 
corresponding to general BA states are "dressed" pseudovacuum eigenvalues, with dressing 
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functions that can be deduced from properties of the transfer matrix such as analyticity, fu- 
sion, etc., independently of the choice of pseudovacuum. (See, e.g., HHj , p2| .) It follows that 



the expressions for the eigenvalues of the transfer matrix corresponding to general BA states 
constructed with pseudovacua um , . . . , um are all the same, and analogously for the pseu- 
dovacua ^a+i) , • • • ,uu\f\- Since the BAE are the conditions that the eigenvalues have van- 
ishing residues, we conclude that the BAE corresponding to the pseudovacua um , . . . , wm 
are all the same; and that the BAE corresponding to the pseudovacua ujn + i) , . . . , un^ are 
all the same. 

We remark that the above proof could be streamlined (in particular, avoiding the argu- 
ment of dressing the pseudovacuum eigenvalues) if one could construct invertible transfor- 
mations U(j) which keep the transfer matrix invariant and change the pseudovacuum, 

t(i) (A , £- , 6+) = % (A , £_ , £+) , 

U {j) U U) = W (J+1) ' ( B - 7 ) 

for j ^ I. For the isotropic case, such transformations are finite elements of the groups SU(l) 
and SU (Af—l), which can easily be constructed 0]. Unfortunately, for the anisotropic case, 



we have not yet succeeded in constructing such transformations for general values of the 
number of sites, iV . Such transformations are presumably related to the so-called g-Weyl 



group 43 , and may have other interesting applications. 



C Symmetry transformations which change the anisotropy 
parameter 

The symmetries of the transfer matrix discussed in Section 2 leave the anisotropy parameter 
rj unchanged. We briefly discuss here related transformations, which however transform 
r) — > —77. Let us now make the dependence on rj explicit, and denote the R matrix ( [2.1 1| ) by 
Ru{^ , T]), and the K matrix ( |2.14|) by Km(X , £ , 77). Evidently, 



R 21 (X,r ] ) = R 12 (X,-r ] ). 
Moreover, let U now be the antidiagonal J\f x M matrix 

( l\ 



(CI) 



u 



(C.2) 



1 



/ 
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One can show that that 



Ut u 2 r 12 (\,v) ur 1 U 2 l 



U K {l) (\,£,r)) U 



-1 



(C.3) 



where V and are given by Eq. ( j2.27| ). The transfer matrix transforms under the corre- 
sponding quantum-space transformation IA as follows: 

u t (l) ( a , f + , n) u- 1 oc t (v) (a , C , e + , -V) • (C4) 

Like the duality transformation (p.26 ), this transformation maps / — > I' and £± — > ^ , but 
it also maps n — > —n. 

Furthermore, define the two M x M matrices 

/ 1 \ / 1 \ 



U, 



(0 



\ 



u, 



(0 



1 / 



1 



V 



(C.5) 



with an antidiagonal upper / x I block, and an antidiagonal lower (A/" — /) x (A/" — /) block, re- 
spectively. The corresponding quantum-space transformations have the following action 
on the transfer matrix 



(C.6) 



However, for A/" > 3, there are not enough of these transformations to determine all the 
elements of the boundary S matrices. 



References 

[1] K. Binder, in Phase Transitions and Critical Phenomena, v. 8, ed. by C. Domb and J. 
Lebowitz (Academic Press, 1983). 

[2] S. Ghoshal and A. B. Zamolodchikov, Int. J. Mod. Phys. A9 (1994) 3841; A9 (1994) 
4353. 

[3] V.V. Bazhanov, Phys. Lett. 159B (1985) 321; Commun. Math. Phys. 113 (1987) 471. 



23 



[4] M. Jimbo, Commun. Math. Phys. 102 (1986) 537; Lecture Notes in Physics, Vol. 246, 
(Springer, 1986), p. 335. 

[5] N.Yu. Reshetikhin, J. Phys. A24 (1991) 2387; CM. Yung and M.T. Batchelor, Nucl. 
Phys. B435 (1995) 430. 

[6] A.E. Arinshtein, V.A. Fateev and A.B. Zamolodchikov, Phys. Lett. B87 (1979) 389; 
A.V. Mikhailov, M.A. Olshanetsky and A.M. Perelomov, Commun. Math. Phys. 79 
(1981) 473; G. Mussardo, Phys. Rep. 218 (1992) 215. 

[7] D. Bernard and A. LeClair, Commun. Math. Phys. 142 (1991) 99. 

[8] V. Pasquier and H. Saleur, Nucl. Phys. B330 (1990) 523. 

[9] P.P. Kulish and N.Yu. Reshetikhin, J. Sov. Math. 23 (1983) 2435; V.G. Drinfel'd, J. Sov. 
Math. 41 (1988) 898; L.D. Faddeev, N. Yu. Reshetikhin and L.A. Takhtajan, Algebraic 
Analysis, 1 (1988) 129; M. Jimbo, Int. J. Mod. Phys. A4 (1989) 3759. 

[10] P.P. Kulish and E.K. Sklyanin, J. Phys. A24 (1991) L435. 

[11] E.K. Sklyanin, J. Phys. A21 (1988) 2375. 

[12] L. Mezincescu and R.I. Nepomechie, J. Phys. A24 (1991) L17. 

[13] L. Mezincescu and R.I. Nepomechie, Mod. Phys. Lett. 40 (1991) 2497. 

[14] L. Mezincescu and R.I. Nepomechie, Int. J. Mod. Phys. 47(1992) 5657. 

[15] H.J. de Vega and A. Gonzalez-Ruiz, Nucl. Phys. 5^7(1994) 553. 

[16] H.J. de Vega and A. Gonzalez-Ruiz, Phys. Lett. B332 (1994) 123, |hep-th/ 9405023 

[17] H.J. de Vega and A. Gonzalez-Ruiz, J. Phys. A26 (1993) L519. 

[18] H.J. de Vega and A. Gonzalez-Ruiz, Mod. Phys. Lett. A9 (1994) 2207. 

[19] M. Jimbo, R. Kedem, T. Kojima, H. Konno and T. Miwa, Nucl. Phys. B44I (1995) 437. 

[20] A. Doikou, L. Mezincescu and R.I. Nepomechie, J. Phys. A31 (1998) 53. 

[21] A. Doikou and R.I. Nepomechie, Nucl. Phys. B521 (1998) 547. 

[22] P.P. Kulish and E.K. Sklyanin, Lecture Notes in Physics, Vol. 151, (Springer, 1982), 
p. 61; V.E. Korepin, N.M. Bogoliubov, and A.G. Izergin, Quantum Inverse Scatter- 
ing Method, Correlation Functions and Algebraic Bethe Ansatz (Cambridge University 
Press, 1993). 



24 



[23] N.Yu. Reshetikhin and M.A. Semenov-Tian-Shansky, Lett. Math. Phys. 19 (1990) 133. 

[24] I.V. Cherednik, Theor. Math. Phys. 61 (1984) 977. 

[25] H.J. de Vega, Int. J. Mod. Phys. A4 (1989) 2371. 

[26] A. A. Belavin, Funct. Anal. Appl. U (1980) 260. 

[27] L.D. Faddeev and L.A. Takhtajan, J. Sov. Math. 24 (1984) 241. 

[28] B. Sutherland, Phys. Rev. B12 (1975) 3795; P.P. Kulish and N.Yu. Reshetikhin, Sov. 
Phys. JETP 53 (1981) 108. 

[29] B.M. McCoy and T.T. Wu, Phys. Lett. 87B (1979) 50; O. Babelon, H.J. de Vega and 
CM. Viallet, Nuc. Phys. B220 (1983) 283. 

[30] M. Grisaru, L. Mezincescu and R.I. Nepomechie, J. Phys. A28 (1995) 1027; A. Doikou, 
L. Mezincescu and R.I. Nepomechie, J. Phys. A30 (1997) L507. 

[31] V.E. Korepin, Theor. Math. Phys. 76 (1980) 165. 

[32] N. Andrei and C. Destri, Nucl. Phys. B231 (1984) 445. 

[33] P. Fendley and H. Saleur, Nucl. Phys. B428 (1994) 681. 

[34] A. Kapustin and S. Skorik, J. Phys. A29 (1996) 1629. 

[35] G. Gasper and M. Rahman, Basic Hypergeometric Series (Cambridge University Press, 
1990). 

[36] A. Forster and M. Karowski, Nucl. Phys. B4O8 (1993) 512; M. Karowski and A. Zapletal, 
J. Phys. A27 (1994) 7419. 

[37] M. Jimbo and T. Miwa, Algebraic Analysis of Solvable Lattice Models, CBMS Regional 
Conference Series in Mathematics vol. 85, (AMS, 1994). 

[38] M. Liischer, Nucl. Phys. Bl 35 (1978) 1. 

[39] A.B. Zamolodchikov, "Fractional-spin integrals of motion in perturbed conformal field 
theory," in Fields, Strings and Quantum Gravity, eds. H. Guo, Z. Qiu and H. Tye, 
(Gordon and Breach, 1989). 

[40] L. Mezincescu and R.I. Nepomechie, Int. J. Mod. Phys. A13 (1998) 2747. 



25 



[41] V.I. Vichirko and N.Yu. Reshetikhin, Theor. Math. Phys. 56 (1983) 805; N.Yu. 
Reshetikhin, Sov. Phys. JETP 57(1983) 691; Lett. Math. Phys. 7(1983) 205; U (1987) 
235. 

[42] L. Mezincescu and R.I. Nepomechie, Nucl. Phys. B372 (1992) 597. 

[43] A.N. Kirillov and N. Reshetikhin, Commun. Math. Phys. 134 (1990) 421. 



26 



